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Abstract. Vertical profiles of the rate of linear change (trend) 
in the altitude range 15-50 km are determined from decadal 
0 3 time series obtained from SCIAMACHY '/EN VIS AT 2 
measurements in limb-viewing geometry. The trends are cal- 
culated by using a multivariate linear regression. Seasonal 
variations, the quasi-biennial oscillation, signatures of the 
solar cycle and the El Nino-Southern Oscillation are ac- 
counted for in the regression. The time range of trend cal- 
culation is August 2002-April 2012. A focus for analy- 
sis are the zonal bands of 20° N-20° S (tropics), 60-50° N, 
and 50-60° S (midlatitudes). In the tropics, positive trends 
of up to 5 % per decade between 20 and 30 km and nega- 
tive trends of up to 10 % per decade between 30 and 38 km 
are identified. Positive O 3 trends of around 5 % per decade 
are found in the upper stratosphere in the tropics and at 
midlatitudes. Comparisons between SCIAMACHY and EOS 
MLS 3 show reasonable agreement both in the tropics and 
at midlatitudes for most altitudes. In the tropics, measure- 
ments from OSIRIS 4 /Odin and SHADOZ 5 are also anal- 
ysed. These yield rates of linear change of O 3 similar to 
those from SCIAMACHY. However, the trends from SCIA- 


1 SCanning Imaging Absorption spectroMeter for Atmospheric 
CHartographY 

2 European environmental research satellite 

3 Earth Observing System (EOS) Microwave Limb Sounder 
(MLS) 

4 Optical Spectrograph and InfraRed Imager System 

5 Southern Hemisphere ADditional OZonesondes 


MACHY near 34 km in the tropics are larger than MLS and 
OSIRIS by a factor of around two. 


1 Introduction 

The stratospheric O 3 layer shields the Earth from UVB and 
UVC radiation in the 240-320 nm range. The absorbed en- 
ergy is the main heat source of the stratosphere, and also 
drives the atmospheric circulation in the upper atmosphere 
(e.g. Weber et al., 2011). The abundance of stratospheric O 3 
is sensitive to both chemical and dynamical factors and has 
been altered by anthropogenic activity. The most prominent 
example of human impact on the O 3 layer is the Antarctic O 3 
hole, which is caused by anthropogenically released chlorine 
and bromine compounds. In order to protect the O 3 layer, the 
production of halogenated source gases was banned by the 
Montreal Protocol of 1987 and its amendments. A slowdown 
of the O 3 decline and even an onset of a recovery have been 
inferred from observations (Newchurch et al., 2003; Yang et 
al., 2006; WMO, 2007). The recovery of the global strato- 
spheric O 3 layer to values observed before the chlorofluo- 
rocarbon (CFC)-driven depletion is, however, not expected 
before the middle of this century (WMO, 2011). 

Measurements of O 3 have a long tradition, and various 
trend studies have been performed using data from the last 
few decades (WMO, 2007). In recent years, several trend 
studies have investigated the impacts of 03 -depleting sub- 
stances (ODS) and dynamical forcings affecting the 0 3 
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shield (WMO, 2007, 2011). The majority of these trend as- 
sessments agreed that the present stage of stratospheric O3 is 
somewhere between the rate of depletion slowing down and 
a subsequent turnaround phase, expected as a result of the 
measures implemented within the Montreal Protocol and its 
amendments. This manuscript contributes to our knowledge 
of stratospheric O3 and its changes during the period follow- 
ing the turnaround in stratospheric halogens in the late 1990s, 
where the halogen load, still high, is slowly declining. 

Vertically resolved O3 trends have been derived from 
observations by satellite instruments, ground-based lidars 
and microwave radiometers, and balloon-borne ozoneson- 
des (Steinbrecht et al., 2006; Jones et al., 2009; Randel 
and Thompson, 2011; Mieruch et al., 2012). In Mieruch 
et al. (2012), O3 trends from SCIAMACHY and those in- 
ferred from a range of instruments were compared within 
the time from 2002 to 2008 with the main focus on possi- 
ble instrumental and algorithmic issues. Analysing a num- 
ber of different satellite datasets, Jones et al. (2009) inferred 
O3 trends from a multi-instrument mean between 1979 and 
2008. SCIAMACHY, OSIRIS/Odin, and SMR/Odin con- 
tributed data starting in the years 2002 and 2001 respec- 
tively. Jones et al. (2009) showed that these satellite data 
agreed with each other to a high extent. After 1997, Jones 
et al. (2009) identified positive midlatitude O3 trends be- 
tween 35 and 45 km on the order of a few percent per decade. 
These findings are corroborated by Steinbrecht et al. (2006, 
2009) by using ground-based measurements. As the upper 
stratosphere is highly sensitive to halogen chemistry (WMO, 
1999), these results may be interpreted as signs of the onset 
of a turnaround in the O3 abundance in response to declin- 
ing halogens. Analysing a merged dataset from satellites and 
ozonesondes, Randel and Thompson (2011) identified nega- 
tive O3 trends in the lowermost tropical stratosphere. These 
were primarily attributed to dynamical factors, namely en- 
hanced tropical upwelling. 

In this study, a global O3 dataset obtained from the SCIA- 
MACHY instrument on-board ENVISAT is investigated. 
SCIAMACHY was in routine operation from August 2002 
until ENVISAT suddenly and unexpectedly failed on 8 April 
2012. In its limb-viewing geometry, SCIAMACHY per- 
formed measurements of the scattered solar radiation. These 
limb measurements are used to retrieve vertical profiles of 
trace gases and aerosols and to detect clouds. 

Rates of linear change or trends in the SCIAMACHY limb 
O3 time series are determined in our investigation for the 
time range from August 2002 to April 2012. These trends 
are presented as functions of altitude and latitude. The result- 
ing vertical trend profiles of O3 are presented for latitudes 
between 60° N and 60° S. The trends from SCIAMACHY 
are compared with those derived from OSIRIS/Odin and 
SHADOZ ozonesondes in the tropics and with those from 
EOS MLS in the tropics and at midlatitudes. 

This paper is structured as follows. Section 2 describes 
briefly the SCIAMACHY limb O3 retrieval. The multivari- 


ate linear regression, used in our study, is introduced and ex- 
plained in Sect. 3. The O3 time series from SCIAMACHY 
and their regression models are described in Sect. 4. In 
Sect. 5, the vertical profiles of the rate of linear change 
or trend of the O3 retrieved from SCIAMACHY measure- 
ments are presented. O3 trend profiles from SCIAMACHY 
are compared with those from EOS MLS and OSIRIS/Odin, 
other atmospheric limb sounders, in Sect. 6. In addition, 
comparisons are performed in the inner tropics including 
SHADOZ ozonesondes. Section 7 discusses atmospheric 
phenomena potentially responsible for the observed O3 
trends. In Sect. 8, sources of errors are discussed. Section 9 
summarises the main results obtained in this study. 


2 SCIAMACHY limb ozone 

The European environmental research satellite ENVISAT 
hosted ten instruments, which were operational until the sud- 
den loss of contact with ENVISAT on 8 April 2012. SCIA- 
MACHY observed electromagnetic radiation upwelling from 
the Earth’s atmosphere in 3 measurement modes: occul- 
tation, nadir, and limb geometry. Detailed information on 
the instrument and mission objectives is provided by Bur- 
rows et al. (1995) and Bovensmann et al. (1999). In limb- 
viewing geometry, the instrument scanned the horizon in 
3.3 km steps from —3 to 92 km (0 to 92 km since October 
2010). This vertical sampling and the instantaneous field of 
view (~ 2.6 km in vertical direction at the tangent point) de- 
fined the vertical resolution, which was typically 3-4 km. 
ENVISAT was in a Sun-synchronous orbit with an incli- 
nation of 98°. This resulted in global coverage for SCIA- 
MACHY limb measurements being achieved within 6 days 
at the Equator and less elsewhere. 

The SCIAMACHY results presented here use the limb O3 
retrieval version 2.9 of IUP Bremen. As an intermediate re- 
trieval result, monthly mean O3 data are gridded horizontally 
into 5° latitude x 15° longitude and vertically into 1 km alti- 
tude. The overall altitude range ofC>3 retrieval is 10 to 75 km. 
As the satellite retrieval often suffered from cloud interfer- 
ence below 15 km, only altitudes above 15 km are used in 
this study. 

Comparing SCIAMACHY limb O3 data to vertically re- 
solved O3 datasets from other instruments, Mieruch et al. 
(2012) showed that comparison results are not significantly 
changed using either collocated measurement pairs or zonal 
monthly means of gridded data. This highlights that zonal 
monthly means are appropriate for performing instrument 
comparisons. Our investigation uses zonal mean O3 data for 
the latitudes defined by the gridding of data as well as for 
wider latitude bands. Determining the rate of linear change 
or trend for all altitudes between 15 and 50 km results in ver- 
tical profiles of O3 trends. 

The zonal monthly mean values are calculated as arith- 
metic means because the errors of single measurements are 
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expected to be normally distributed and no issue with out- 
liers is known. Any of these zonal monthly means is typi- 
cally composed of several hundreds of single measurements. 
Hence, the random errors of zonal monthly means can be 
deemed as negligible. Trend uncertainties used in the follow- 
ing are calculated from the covariance matrix of regression. 

The boundaries of 60° N and 60° S are chosen for this in- 
vestigation to avoid any effects which are directly related to 
the polar vortex and because of gaps in the SCIAMACHY 
sampling during the high-latitude winter. In order to avoid 
possible influences of the South Atlantic Anomaly, measure- 
ments with tangent point coordinates between longitudes of 
90° W and 0° W and latitudes south of 20° S are excluded 
from this analysis (it is worth noting that a typical difference 
in latitude between instrument and tangent point coordinates 
of around 20 ° is taken into account by the criteria used). 

3 Trend determination 

3.1 Multivariate linear regression 

A multivariate linear regression is applied to the datasets with 
the autocorrelation of consecutive values accounted for by 
the Cochrane-Orcutt transformation (Cochrane and Orcutt, 
1949). The observed values of autocorrelation are similar to 
those shown in Mieruch et al. (2012). The regression is un- 
weighted. The following terms are used in the regression: 

fj. + a)t + J2 ^$l;Sm(^^) + $2 7 cos(^^)^ + R(Q, (1) 

where t is the time in months and n, a>, and 
$ 1 i , ySi 2 , . .., $23, $24 are the fitting parameters. The 
harmonics with 12-, 6 -, 4-, and 3-month periods, j = 1, 2, 3, 
4, a sine and cosine term each, are used to represent seasonal 
changes. The harmonics of 12 and 6 months approximate 
the annual and semi-annual cycles. The combination of sine 
and cosine terms provides a flexible adjustment to any phase 
of the (semi-)annual variation. As discussed by Stiller et al. 
(2012), the terms having 4- and 3-month periods improve 
the fit quality because the oscillation patterns are possibly 
compressed or stretched with respect to a harmonic shape. 
In addition, a linear fit having a gradient or slope a> and 
an intercept /r is included in the regression analysis. The 
gradient represents the O 3 trend at a given altitude. The 

1 <r value defined by the covariance matrix of regression 
coefficients serves as trend uncertainty. It is also used to 
draw error bars of trends in subsequent figures. The trends 
and their uncertainties are presented in relative units of 
% dec -1 , i.e. percent per decade (relative with respect to the 
mean value of the underlying time series). The criterion for 
the trend being significant at the 95 % confidence level is that 
the absolute ratio of the trend to its uncertainty is larger than 

2 (Tiao et al., 1990). As discussed later in Sect. 3.2, it is also 


advantageous to extend the linear regression (see Eq. 1) by 
additional terms representing the quasi-biennial oscillation 
(QBO), the solar cycle (SC), and the El Nino-Southern 
Oscillation (ENSO): 

R(t) = QBO(f) + SC(f) + ENSO(f). (2) 

At northern midlatitudes, Eq. (1) is optimised because of the 
strong interannual variability of the annual cycle of O 3 . In 
the 15-26 km range of the 60-50° N latitude band, the cu- 
mulative eddy heat flux (based on ERA-lnterim, 50 hPa eddy 
heat flux integrated from 45° N to 75° N, cumulation starts 
in October) is fitted instead of harmonic terms. The cumula- 
tive eddy heat flux is a proxy accounting for the wave-driven 
O 3 accumulation in winter and early spring. For 60-50° N, 
the interannual variability in total O 3 is found to be well ap- 
proximated by the cumulative eddy heat flux (Dhomse et al., 
2006). Here, correlation tests show that throughout the al- 
titude range of 15-26 km, the proxy being time lagged by 
2 months is optimum. 

3.2 QBO, solar cycle, and ENSO signatures in ozone 

The quasi-biennial oscillation (QBO) is a quasi-periodic sig- 
nal observed in the tropical stratospheric zonal wind speed. 
It alternates between east and west phase with the period 
varying between 2 and 3 yr. Previous studies report a max- 
imum QBO response of tropical O 3 in the altitude ranges of 
20-27km and 30-38km (Zawodny and McCormick, 1991; 
Chipperfield et al., 1994). 

In the tropics, the QBO response of O 3 is driven by ef- 
fects related to the vertical transport (Butchart et al., 2003). 
The QBO-induced meridional circulation modulates the up- 
welling branch of the Brewer-Dobson circulation, leading to 
the largest changes of the O 3 volume mixing ratio (VMR) in 
the altitude range of large vertical gradients of the O 3 VMR, 
which are located below and above the maximum in the ver- 
tical O 3 VMR profiles. The result is a compacting (vertical 
squeezing) of the mean vertical O 3 VMR profile at west- 
erly QBO phase and a lofting (vertical stretching) at east- 
erly QBO phase. Around the maximum of the mean vertical 
profile of the O 3 VMR, where the vertical gradient changes 
sign, the phase change of the QBO response of O 3 takes place 
(Butchart et al., 2003). With the QBO phases modulating the 
Brewer-Dobson circulation, a QBO response is expected to 
be observed also in extratropical O 3 (Baldwin et al., 2001). 

In the following, the lOhPa and 30hPa Sin- 
gapore winds (Naujokat, 1986, available from 
http://www.geo.fu-berlin.de/en/met/ag/strat/produkte/qbo/) 
are used as fit proxies: 

QBO(f) = aQBOio(f) + bQBO 30 (t), (3) 

with a and b being their fit coefficients. The indices “10” and 
“30” denote the pressure level of the proxy. This approach 
is applied both in the tropics and extratropics to account for 
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any QBO effects. The monthly time series of the Singapore 
winds are smoothed by a 4-month running average before be- 
ing used as QBO proxy. The same procedure is done with the 
Mg II index when used as a solar cycle proxy. Both proxies 
are not detrended. 

The 1 1 yr solar cycle, which leads to pronounced varia- 
tions of the solar radiation in the U V spectral range, is known 
to influence stratospheric O 3 (Gray et al., 2010). It is widely 
accepted that the solar cycle response of O 3 occurs without 
any time lag and at positive correlation (e.g. Soukharev and 
Hood, 2006; Remsberg and Lingenfelser, 2010). The solar 
cycle response of O 3 may extend from the tropics into the 
midlatitudes (SPARC CCMVal, 2010). In the following, the 
solar cycle proxy is included in the multivariate linear regres- 
sion for any latitude analysed. A multi-instrument monthly 
average Mg II index - which is mainly based on GOME, 
SCIAMACHY, and GOME-2 observations over the time 
span under study - serves as solar cycle proxy. The under- 
lying data can be obtained from http://www.iup.uni-bremen. 
de/gome/solar/MgII_composite.dat (Snow et al., 2013; We- 
ber et al., 2013). 

In addition to the QBO and the solar cycle, the O 3 time 
series may also be influenced by the El Nino Southern Oscil- 
lation (ENSO) (Randel et al., 2009; Randel and Thompson, 
2011). Employing an ENSO proxy based on the Nino 3.4 
index, Oman et al. (2013) previously demonstrated its ap- 
plicability to reveal ENSO signatures in vertically resolved 
O 3 time series from EOS MLS. With the ENSO response of 
O 3 propagating towards higher altitudes, time lags between 
ENSO signatures in O 3 and ENSO proxies may arise. For 
instance, Randel et al. (2009) and Randel and Thompson 
(201 1) employed ENSO proxies at a time lag of 2 months. 

Similar to Oman et al. (2013), an ENSO proxy based on 
anomalies of the Nino 3.4 index (available from http://www. 
cpc.ncep.noaa.gov/data/indices/) is used in our study. Time 
lags are accounted for by the combined fit of the proxy and 
its time derivative: 


ENSO(f) = c 


/V34(0 + 


d/V 34 (0 
~ dr 


A(r) , 


(4) 


where /V 34 denotes the proxy and A (t ) the time lag. Account- 
ing for the non-linearity of the fit is repeated with the 

lagged proxy until the remaining time lag approaches zero. 
This iteration starts at a time lag of 2 months. The time lag of 
the proxy not settling at a certain value (including zero) or be- 
coming unrealistic is the criterion for omitting ENSO terms 
at regression. In practice, the latter is the case at extratropical 
latitudes and in the tropics at altitudes above 25 km. In the 
following, ENSO is considered at regression only at tropi- 
cal latitudes (any latitude band within 20° N-20° S) in the 
15-25 km altitude range. If the criterion for omitting ENSO 
terms is still fulfilled, the underlying time series is excluded 
from analysis (resulting in the gaps in the vertical trend pro- 
files shown Figs. 8 and 11 and in the QBO amplitude/solar 
cycle response shown in Figs. Al and A2). 


4 Time series of ozone 

Figures 1-6 show O 3 time series at the altitudes of 21 km, 
35 km, and 44 km at northern midlatitudes (60-50° N) and 
in the tropics (20° N-20° S). The time series are overlaid 
by their fitting curves from the regression model, and the 
fit residuals are plotted below. The individual terms from 
the regression - i.e. the harmonic, QBO, solar cycle, and 
linear terms - are shown in corresponding panels below 
(in the tropics at 21 km, ENSO terms are included as well 
while, at northern midlatitudes at 2 1 km, eddy heat flux terms 
are shown instead of harmonic terms). The harmonic terms 
are shown comprising all periods used, i.e. 12, 6 , 4, and 3 
months, and also as combination of the terms with 4- and 
3-month periods only. With the exception of 44 km at 60- 
50° N, the latter terms of 4- and 3-month periods improve 
the closure between time series and fitting curve, but do not 
change the amplitude of the harmonic terms substantially. 
The harmonic terms are overlaid by the time series with all 
fit terms with the exception of the harmonics subtracted. Any 
other component of the fitting curve is overlaid by the time 
series with all other fit terms removed accordingly. 

At all altitudes, good fit quality is obtained. The use of 
harmonics for approximating the (semi-)annual oscillation as 
well as the combination of the lOhPa and 30hPa Singapore 
winds to fit QBO appears to be a sufficiently complete ap- 
proach. The fit residuals are well within ±5% to ±10% of 
the fitted time series both at midlatitudes and in the tropics. 

At 21 and 35 km at midlatitudes, the annual cycle terms 
(green curve has an annual pattern) dominate in the fitting 
curve. At 21km, O 3 has its annual maximum in late win- 
ter/early spring and declines throughout summer. At 35 km, 
the annual maximum of O 3 is reached in late spring/early 
summer, persisting for the summer months. At 44 km at mid- 
latitudes, the harmonic terms have two maxima and minima 
per year. The annual minimum occurs in summer and a sec- 
ondary minimum in winter. These patterns are explained as 
follows: at 21 km, the periodic variability of O 3 is dominated 
by transport effects in winter, leading to an accumulation of 
O 3 , and by catalytic O 3 depletion in summer which is re- 
sponsible for the descending branch of the annual cycle. At 
35 km and 44 km, the periodic variability is governed by pho- 
tochemical O 3 production and by catalytic O 3 depletion re- 
spectively, both being at its maximum in summer (Perliski et 
al., 1989). At 35 km, the annual maximum in summer driven 
by photochemical production is flattened by catalytic deple- 
tion which is of minor consequence. At 44 km, catalytic de- 
pletion has gained control leading to the annual minimum in 
summer. Overlaid photochemical production of O 3 leads to a 
secondary minimum in winter. 

At 21km in the tropics, the annual cycle and the QBO 
make comparable contributions to the fitting curve. Dis- 
cernible contributions are also made by the ENSO. The an- 
nual cycle has its maximum in the Northern Hemisphere 
summer. It resembles well the phase of the annual oscillation 
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Fig. 4. Same as Fig. 1 but for 35 km altitude and 20° N-20° S. Third panel as in Fig. 2: harmonics as seasonal terms. 


of tropical O 3 in the lowermost tropical stratosphere which is 
driven by effects related to the tropical upwelling (Randel et 
al., 2007). At 35 km in the tropics, the QBO and semi-annual 
terms dominate the periodic variability of O 3 . As discussed 
before, the maximum QBO response of tropical O 3 is located 
within around 20-27 km and 30-38 km. At 44 km in the trop- 
ics, the semi-annual terms provide the most important contri- 


bution to the fitting curve. In the middle tropical stratosphere, 
the semi-annual oscillation of O 3 is driven by photochemical 
O 3 production, which reflects two Equator passages of the 
Sun per year (Perliski et al., 1989). 
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Fig. 5. Same as Fig. 1 but at 44 km altitude. Third panel as in Fig. 2: harmonics as seasonal terms. 
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Fig. 6. Same as Fig. 1 but for 44 km altitude and 20° N-20° S. Third panel as in Fig. 2: harmonics as seasonal terms. 


5 Ozone trends from SCIAMACHY 

O 3 trend profiles and estimated uncertainties (both in 
% dec -1 ) are shown in Fig. 7 for the period of August 2002- 
April 2012. The underlying latitude bands are 60-50° N and 
50-60° S (midlatitudes) and 20° N-20° S (tropics). The fol- 
lowing features of particular interest are readily identified: 


- double-peak structure of positive trend in the tropics 
below 30 km (having two maxima at values close to 
5 % dec - 1 between 20 and 30 km) 

- single-peak structure of negative trend in the tropics 
between 30 and 38 km (maximising at 33-34 km at a 
value of around — 10 % dec -1 ) 
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Fig. 7. Ozone trend profiles inferred from SCIAMACHY limb 
measurements for the latitude bands 60-50° N, 20° N-20° S, and 
50-60° S, August 2002-April 2012. Error bars: Iff. 

- positive trends in the upper tropical stratosphere be- 
tween 39 and 47 km and at northern midlatitudes (be- 
tween 40 and 45 km) of around 5 % dec -1 

- hemispheric asymmetry of trends at midlatitudes be- 
tween 25 and 35 km: negative trends in the north and 
zero to positive trends in the south maximising at val- 
ues of around — 10%dec -1 and 3%dec -1 respec- 
tively. 

In the 35-45 km range of 50-60° S, the trends are rela- 
tively close to zero and become increasingly negative above 
45 km. The trend errors are comparable for all three latitude 
bands. 

A stability check was performed for the O3 trend profiles 
from SCIAMACHY by shifting the starting point of the un- 
derlying time series forward in time in 6-month steps, i.e. 
from August 2002 to February 2003, from February 2003 to 
August 2003, from August 2003 to February 2004, and from 
February 2004 to August 2004. In all cases, the trends and 
their uncertainties are well replicated and the overall shape 
of vertical trend profiles is conserved for the three latitude 
bands. It is, however, essential to extend the period to 2011. 
Otherwise, the solar cycle terms reflect mainly the descend- 
ing branch of the solar cycle and may interfere with the lin- 
ear trend. This may, in turn, result in larger trend uncertain- 
ties. Practically no change of the vertical trend profiles from 
SCIAMACHY results if omitting data from the year 2012, 
i.e. confining the time range of trend calculation to August 
2002-December 2011. 

Figure 8 shows a latitude-altitude plot of the O3 trends 
spanning the latitude range from 60° N to 60° S. The plot is 
based on zonally averaged O3 data in 5° latitude steps, de- 
fined by the data gridding described in Sect. 2. As a function 
of the latitude, the trends are shown for the 15 to 50 km alti- 
tude range. As evidenced before for the 20° N-20° S tropical 


%/dec 
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Fig. 8. Latitude-altitude plot of ozone trends from SCIAMACHY: 
latitude range from 60° N to 60° S, altitude range from 15 to 50 km. 
Non-hatched are significances at the 2cr level. August 2002-April 
2012. Legend on the right: units of % dec -1 . 

latitude band, there are double-peaked positive trends below 
30 km (significant at their maxima at around 21 km and be- 
tween 25 and 30 km) and single-peaked negative trends be- 
tween 30 and 38 km in the inner tropics. In the 30-35 km 
range, an “island” of significant negative O3 trends of up to 
—20 % dec -1 now appears. At latitudes from 20° N to 30° N, 
significant positive trends are seen at altitudes between 20 
and 25 km. Positive trends in the upper stratosphere (iden- 
tified in Fig. 7 for 20°N-20° S between 39 and 47 km and 
also for 60-50° N between 40 and 45 km) are observed over 
large parts of the latitude range from 60° N to 60° S. These 
trends are mostly significant at values of around 5% dec -1 . 
The hemispheric asymmetry at midlatitudes between 25 and 
35 km is clearly pronounced with significant negative trends 
in the Northern Hemisphere and significant positive trends in 
the Southern Hemisphere. At around 15 km, positive trends 
appear between 20° S and 30° S. These trends should, how- 
ever, be treated with caution because of potential impact of 
clouds at altitudes below 20 km (Mieruch et ah, 2012). 

A discussion of the latitude and altitude dependence of the 
QBO and solar cycle responses as seen by SCIAMACHY 
limb is presented in Appendix A. 

6 Comparisons of trends from SCIAMACHY with 
those from other instruments 

6.1 OSIRIS/Odin 

The Odin satellite, which carries the OSIRIS instrument 
(Llewellyn et ah, 2004), follows a near-terminator orbit so 
that the illumination conditions specific to OSIRIS are al- 
ways close to the solar terminator. The ascending node was 
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between 6 and 7 p.m. over the time span considered here. 
Similar to SCIAMACHY, OSIRIS measures the scattered so- 
lar light in limb-viewing geometry. Seasonal variations of 
day- and nighttime interrupt the OSIRIS measurements in 
the winter hemisphere. As a result, OSIRIS provides con- 
tinuous time series only in the tropics. The data might also 
be inhomogeneously distributed around the Equator, i.e. pre- 
dominantly to the north from it during northern summer and 
vice versa. Thus, the comparison between the instruments is 
performed for the tropical 20° N-20° S latitude band only. 

The OSIRIS O 3 data, used in this study, are from 
the retrieval version 5.07 produced by the University of 
Saskatchewan, which has a modified filtering of outliers with 
respect to earlier versions. The OSIRIS data separate clearly 
into records at a.m. and p.m. local time. As pointed out by 
Huang et al. (2010), diurnal variations of O 3 are expected 
to be on the order of only a few percent in the stratosphere. 
The AM Equator crossing time of Odin is closer to the Equa- 
tor crossing of ENV1SAT at daylight. Accordingly, the a.m. 
data from OSIRIS, which are closer in local time to the data 
from SCIAMACHY, are used in the following (without any 
correction for diurnal variation). 

The retrieval altitudes of OSIRIS and SCIAMACHY lie 
between one another. For each altitude of SCIAMACHY 
from 15 to 50 km, the O 3 values from the adjacent altitudes 
of OSIRIS were interpolated. 

The instrument comparison between SCIAMACHY and 
OSIRIS (also between SCIAMACHY and EOS MLS as de- 
scribed in Sect. 6.2) uses zonal monthly mean data. As men- 
tioned in Sect. 2, a comparison of zonal monthly means from 
gridded data is considered an equivalent alternative to the 
use of collocated data. As a further check, the trend profiles 
from SCIAMACHY and OSIRIS (comparison described in 
Sect. 6.4) were tested to remain practically unchanged when 
using collocated measurements. 

6.2 EOS MLS 

The MLS instrument aboard the EOS Aura satellite (Jarnot 
et al., 2006) measures the thermal radiation emitted by the 
Earth’s atmosphere in the microwave spectral range on both 
the day and night side of the near-polar orbit. MLS went 
into routine operation in August 2004. In this comparison, 
we have used O 3 profiles retrieved from MLS daytime mea- 
surements. Aura crosses the Equator in ascending node at 
13:45 local time. Its daytime measurements extend deep into 
both hemispheres similar to those from SCIAMACHY. Thus, 
the comparison between the instruments is also done at mid- 
latitudes, represented by the latitude bands of 60-50° N and 
50-60° S. The period used for the comparison is August 
2004-April 2012. 

In accordance with the recommendations from the MLS 
team, the MLS O 3 data used here are from the retrieval ver- 
sion 2.2. The pressure serves as the vertical coordinate in ver- 
sion 2.2. After conversion from pressure to geometric height, 


the MLS O 3 profiles are interpolated at the altitudes of the 
SCIAMACHY retrieval. Zonal monthly means are calculated 
between 15 and 50 km. 

6.3 SHADOZ ozonesondes 

Balloon-borne in situ measurements of ozonesondes are also 
used for a comparison with the O 3 trend from SCIAMACHY 
limb. The altitude range of the balloon sonde O 3 data is lim- 
ited to 30-35 km by different factors, such as balloon burst 
and evaporation and/or freezing of the wet-chemical sensing 
solution used in the sondes. The data from ozonesondes have 
a vertical resolution on the order of 100 m, clearly finer than 
that from satellite instruments. For the comparison, the ver- 
tical O 3 profiles from the sondes are smoothed over intervals 
of ±2 km relative to altitudes defined by the O 3 retrieval from 
SCIAMACHY, i.e. from 15 to 30 km in steps of 1 km. In the 
following, these smoothed O 3 profiles are used to determine 
monthly averaged O 3 time series from the sondes. 

SHADOZ is a network of more than 10 ozonesonde sta- 
tions at (sub)tropical sites (e.g. Thompson, 2003). Since 
coming into being in 1998, the network has been systemati- 
cally operational with typically several ozonesonde launches 
per month per station. The datasets of all stations are 
freely accessible via http://croc.gsfc.nasa.gov/shadoz/. Here, 
we select the inner-tropical stations at Ascension (8.0° S, 
14.4° W), Kuala Lumpur (2.7° N, 101.7° E), Nairobi (1.3° S, 
36.8° E), Natal (5.4° S, 35.4° W), and Paramaribo (5.8° N, 
55.2° W). These stations are relatively equidistant in latitude 
(less importance is attached to their longitudes as tropical 
stratospheric O 3 is known to be highly zonally symmetric). 
The monthly averaged O 3 time series over all 5 stations are 
calculated. In order to account for possible offsets between 
the stations, the mean values of single-station time series are 
forced to be equal before averaging. The is done by shift- 
ing them so that a reference value is matched. This reference 
value is the mean value of one of the stations, in our case 
Paramaribo (arbitrarily selected). The resulting sonde trend 
profile is compared to the trend profile from SCIAMACHY, 
which is based on data selected at latitudes and longitudes 
around the sites of the ozonesonde stations and averaged ac- 
cordingly. 

6.4 Comparison results 

SCIAMACHY vs. OSIRIS and MLS 

Figures 9 and 10 show the comparison between SCIA- 
MACHY and MLS at midlatitudes for the period August 
2004-April 2012. At northern midlatitudes, the MLS O 3 
trend profile has less variability with altitude than that of 
SCIAMACHY. There are some significant deviations be- 
tween 30 and 35 km, with SCIAMACHY showing nega- 
tive trends while MLS shows small positive trends. Positive 
trends seen by SCIAMACHY around 40 to 45 km agree with 
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Fig. 9. Comparison of ozone trends inferred from SCIAMACHY 
and MLS measurements: 60-50° N, August 2004-April 2012. Error 
bars: lcr. 

MLS. There is also reasonable agreement at southern mid- 
latitudes. The trend profiles from both instruments follow 
each other closely and have similar error bars. The agreement 
between SCIAMACHY and MLS holds up to 50 km, where 
both profiles become increasingly negative. 

Tropical 0 3 trends from SCIAMACHY, OSIRIS, and 
MLS are compared in Fig. 11. The tropics are represented 
by the latitude band of 20° N-20° S, and the time span from 
August 2004 to April 2012 is selected for the comparison. 
In qualitative agreement, SCIAMACHY, OSIRIS, and MLS 
show negative trends between 30 and 35 km. SCIAMACHY 
reaches around — 10%dec _1 , while the trends from OSIRIS 
are around — 4%dec _1 at maximum. The trends from MLS 
lie between those from SCIAMACHY and OSIRIS. 

Between 25 and 30 km, the O3 datasets from all three in- 
struments show in good agreement positive trends. Agree- 
ment within the trend uncertainties is achieved down to 
around 20 km. In the UTLS region, the trends from SCIA- 
MACHY and OSIRIS change from moderately positive at 
20 km to moderately negative at 15 km, while MLS increases 
to 10-20% dec -1 (which is consistent with the results from 
SHADOZ; see following subsection). Between around 38 
and 46 km, both SCIAMACHY and OSIRIS show positive 
trends. Trends up to around 6% dec -1 are obtained from 
SCIAMACHY near 45-46 km and from OSIRIS around 
40 km. Contrastingly, the trends from MLS are zero to nega- 
tive above 40 km. Above 48 km, zero to slightly negative O3 
trends are shown by all compared instruments. Taken as a 
whole, the shape of the trend profile of tropical O3 identified 
from SCIAMACHY in Sect. 5 is confirmed by both OSIRIS 
and MLS. 



trend in %/dec 

Fig. 10. Comparison of ozone trends inferred from SCIAMACHY 
and MLS measurements: 50-60° S, August 2004-April 2012. Error 
bars: lcr. 

SCIAMACHY vs. SHADOZ 

Figure 12 compares the SHADOZ sonde and SCIAMACHY 
O3 trend profiles from 15 to 30 km. Though it is derived 
from data averaged over latitudes and longitudes matching 
the sites of ozonesonde stations, the trend profile from SCIA- 
MACHY resembles well the trend pattern derived from zonal 
averages at latitudes from around 5° N to 5° S in Fig. 8. This 
is a result of the high zonal symmetry of tropical O3 (Thomp- 
son, 2003; Randel and Thompson, 2011). The trend profile 
from SCIAMACHY shows zero to positive trends between 
15 and 30 km. Good agreement to the positive trends seen by 
the sondes is obtained between 18 and 29 km. Below 18 km, 
the trends from the sondes climb towards a maximum while 
the trends from SCIAMACHY decay towards zero. The dif- 
ferences between 15 and 20 km may arise from the satel- 
lite retrieval suffering from the impacts of high convective 
clouds; see e.g. Mieruch et al. (2012). 

7 Discussion of factors potentially explaining the O 3 
trends observed 

The O3 trend profiles from SCIAMACHY and the compari- 
son instruments tend to have non-negative trends in the trop- 
ical 17-20 km range. Constraining our study of trends to the 
2002 to 2012 decade, this behaviour is in contrast to the re- 
sults obtained by Randel and Thompson (2011). Their O3 
trends were inferred from a dataset merged between SAGE II 
and SHADOZ spanning the time range of 1984-2009. Neg- 
ative trends were identified in the lowermost tropical strato- 
sphere and attributed to enhanced tropical upwelling. 

In the tropical stratosphere below 30 km, both changing 
dynamics and chemical composition potentially explain the 
positive O3 trends observed by SCIAMACHY For exam- 
ple, changes in the rate of tropical upwelling, changes in 
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Fig. 11. Comparison of ozone trends inferred from SCIA- 
MACHY, MLS, and OSIRIS measurements: 20°N-20°S, Au- 
gust 2004-April 2012. Error bars: lcr. 


the amount of halogens, or changes in the NOx amount are 
possibly responsible for changes in O 3 . Analysing O 3 par- 
tial columns between 18 and 25 km, Yang et al. (2006) found 
that halogens are of minor importance in this altitude range 
because CFCs do not release radicals immediately after their 
injection through the tropical tropopause. Since the vertical 
gradient of O 3 is positive, enhanced vertical advection has a 
negative O 3 response in the lower tropical stratosphere. Ran- 
del and Thompson (2011) suggested that an increase of the 
tropical upwelling has led to negative O 3 trends around 17 to 
21 km observed over the period 1984-2009. As the O 3 trends 
seen by SCIAMACHY and the comparison instruments are 
non-negative at this altitude range and become clearly pos- 
itive between 25 and 30 km, they must have another expla- 
nation. According to Nevison et al. (1999), an increase in 
NOx might be a reason for a positive trend in O 3 at these 
altitudes. This is because NOx is involved in buffering ef- 
fects as a result of interactions with HOx and ClOx beside 
catalytic O 3 depletion. Nevison et al. (1999) have shown that 
an increase in NOx results in an O 3 enhancement extend- 
ing up to around 30 km in the proximity of the Equator. As 
noted by Burrows et al. (1979), NOx may also catalytically 
produce O 3 in the lowermost stratosphere. The altitude range 
below 30 km coincides with the positive O 3 trends observed 
in our study. As also follows from Nevison et al. (1999), in- 
creasing NOx may lead to O 3 decreases above 30 km. Ac- 
cordingly, increasing NOx may contribute to the negative O 3 
trends seen by SCIAMACHY in the tropical 30-38 km range 
as well. However, the anthropogenic increase of NOx as- 
sumed in Nevison et al. (1999) results in a decrease of around 
1 % per decade in tropical midstratospheric O 3 , which cannot 
explain the observed negative trends of —4 to — 10%dec _1 
inferred in our study (cf. Fig. 11). Thus, the correlation be- 
tween O 3 and NOx pointed out by Nevison et al. (1999) can 
explain the observed trends only if sources of NOx other than 


Fig. 12. Comparison of ozone trends inferred from SCIAMACHY 
and from an average of tropical SHADOZ ozonesonde stations 
including Ascension (8.0° S, 14.4° W), Kuala Lumpur (2.7° N, 
101.7° E), Nairobi (1.3° S, 36.8° E), Natal (5.4° S, 35.4° W), and 
Paramaribo (5.8° N, 55.2° W). The data from SCIAMACHY are se- 
lected at latitudes and longitudes close to the sites of the ozonesonde 
stations and averaged accordingly. The considered period is August 
2002-April 2012. Error bars: la. 


anthropogenic emissions exist. An investigation of the possi- 
ble sources of NOx is beyond the scope of this study. Further 
research is needed on this issue. 

Positive O 3 trends of approximately 5% dec -1 are in- 
ferred from SCIAMACHY in the upper stratosphere over 
large parts of the tropics and midlatitudes. As O 3 in the upper 
stratosphere is highly sensitive to halogens (WMO, 1999), 
these trends are favoured to reflect declining halogens. Ev- 
idence for a turnaround in upper stratospheric O 3 has al- 
ready been pointed out by Reinsel (2002). Some attribution 
to declining halogens has been provided by Newchurch et 
al. (2003). Their study is based on satellite data indicating 
at least a levelling-off of O 3 in the tropics and extratropics. 
Meanwhile the O 3 abundance continuing some levelling-off 
or slightly increasing is supported by various ground-based 
lidar and microwave radiometer stations (Steinbrecht et al., 
2006, 2009). 

A potential explanation for some of the O 3 trends obtained 
as a function of latitude and altitude is provided by Eckert 
et al. (2013) based on MIPAS/ENVISAT O 3 data covering 
the same observation period as SCIAMACHY. A hypotheti- 
cal shift of subtropical mixing barriers is theorised to result 
in negative trends in the northern extratropics and positive 
trends in the southern extratropics. This pattern of trends is 
similar to the asymmetry between the northern and southern 
trend profiles from SCIAMACHY between 25 and 35 km. 
In addition, positive/negative O 3 trends in the lower/middle 
tropical stratosphere are theorised by Eckert et al. (2013), 
although the tropical trends from SCIAMACHY partly are 
more pronounced and somewhat different in altitude. This 
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explanation, if correct, implies that the subtropical mixing 
barriers having undergone shifts during the last decade may 
contribute to the tropical O 3 trends observed in our study. 


8 Potential instrumental and algorithmic issues 

As discussed by von Savigny et al. (2005), pointing misalign- 
ment is one of the major error sources for limb-viewing in- 
struments on ENV1SAT. A robust method for determining 
tangent height offsets from solar occultation measurements 
of SCIAMACHY has been developed by Bramstedt et al. 
(2012). Using this method, tangent height offsets on the order 
of few hundred metres were identified. The overall behaviour 
of the pointing misalignment is mainly governed by constant 
and seasonally varying components. In addition, a continu- 
ous drift of the pointing misalignment on the order of 50 m 
throughout the first four years of the SCIAMACHY mission 
and no drift thereafter were identified. While any possible 
artefacts due to constant and seasonally varying components 
are neutralised by the corresponding regression terms and 
hence are not expected to produce trend errors, the influence 
of the continuous drift on the resulting trends needs to be fur- 
ther investigated. Assuming similar pointing drifts for both 
limb and occultation measurements, this is done by calcu- 
lating the trends for drift-corrected time series. The weight- 
ing functions of the vertical profile of O 3 with respect to the 
pointing offset needed to perform the correction are calcu- 
lated numerically using one representative profile for each 
season. The corrected time series yield up to 1 % dec -1 lower 
O 3 trends in the 40-50 km range and practically the same re- 
sults as uncorrected elsewhere. 

As described in Sect. 6.4, negative O 3 trends from SCIA- 
MACHY observed between 30 and 38km in the trop- 
ics are qualitatively confirmed by both EOS MLS and 
OSIRIS/Odin, but the comparison instruments show less 
negative trends than SCIAMACHY. This altitude range is 
known to be possibly affected by reduced sensitivity in the 
SCIAMACHY retrieval resulting from a transition between 
the retrieval information coming from UV and visible wave- 
lengths. Previously, positive offsets in SCIAMACHY O 3 
number densities compared to the datasets from other instru- 
ments were observed in the tropics around 30 km and around 
40 km by Mieruch et al. (2012). Similarly, compared to a va- 
riety of satellite instruments in Tegtmeier et al. (2013), SCIA- 
MACHY has 10-20 % higher O 3 mixing ratios in the tropics 
around lOhPa, which is close to 30 km in altitude. In addi- 
tion, the QBO amplitude seen by SCIAMACHY lies above 
the multi-instrument average in the middle tropical strato- 
sphere. However, the QBO amplitude, if overestimated, is not 
necessarily a source of error for trends as QBO signals influ- 
ence mainly trend uncertainties. As a result of the inclusion 
of QBO terms into the regression model used in this study, 
trend uncertainties of SCIAMACHY are clearly reduced in 
the middle and lower tropical stratosphere while trend values 


remain practically unchanged. Furthermore, analysing long- 
term time series, Tegtmeier et al. (2013) found no indications 
for temporal instability of the SCIAMACHY instrument. 

As pointed out by Jones et al. (2009) and McLinden and 
Fioletov (2011), O 3 trends can differ dependent on the phys- 
ical units used for ozone and the vertical coordinate (VMR 
versus number density, pressure versus geometric height) due 
to possible trends in the neutral density (temperature and 
pressure) that is needed for unit conversion. Our investigation 
shows that the O 3 trends in number density units, obtained 
from EOS MLS measurements using their own pressure and 
temperature information, agree better with SCIAMACHY in 
the tropical 30-40 km range, while no significant changes 
are observed below 30 km (not shown). However, converting 
number density into volume mixing ratio does not alter the 
ozone trends for OSIRIS/Odin and SCIAMACHY. This sug- 
gests that the neutral density trends derived from ECMWF 
data and those obtained from MLS are significantly different 
at these altitudes. Both SCIAMACHY and OSIRIS do not 
retrieve temperature (unlike MLS), and, therefore, ECMWF 
operational temperatures are used in the retrieval and unit 
conversion. As pointed out by Adams et al. (2013), drifts in 
the meteorological data can also affect ozone trends when 
units are converted. In our study, a significant disagreement 
above 35 km was identified when comparing the neutral den- 
sity trends obtained from ECMWF operational and ERA- 
Interim datasets (not shown here). 

9 Conclusions 

Rates of linear change or trends were determined for strato- 
spheric O 3 in the 15-50 km range for the time range Au- 
gust 2002 to April 2012. In the tropics (20°N-20°S), the 
O 3 trends from SCIAMACHY are found to be positive in the 
lower stratosphere and negative in the middle stratosphere. 
The positive trends maximise between 20 and 30 km, are 
double-peaked, and maximise at values close to 5% dec -1 . 
The negative trends are single-peaked and maximise at 33- 
34 km at values of around — 10 % dec - 1 . The observed trends 
at this altitude are about twice the trends observed by OSIRIS 
(this paper), MIPAS (Eckert et al., 2013), and the combined 
SAGE II/GOMOS (Kyrola et al., 2013). 

As suggested by Nevison et al. (1999), increasing NOx 
may contribute both to the positive O 3 trends observed be- 
low 30 km and the negative O 3 trends observed above 30 km 
in the tropics. In the upper stratosphere between around 40 
and 47 km, positive O 3 trends of around 5% dec -1 are seen 
by SCIAMACHY from the tropics to the midlatitudes of both 
hemispheres. These trends are believed to be in response to 
declining halogens as this region is highly sensitive to halo- 
gen chemistry (WMO, 1999). Similar O 3 trends have been 
previously reported by Steinbrecht et al. (2006, 2009), Jones 
et al. (2009), and WMO (2011) in the 35-45 km range. 
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For the latitude band of 20° N-20° S, representing the 
tropics, O 3 trends have been compared between SCIA- 
MACHY, OSIRIS/Odin, and EOS MLS. In addition, a 
trend comparison between SCIAMACHY and SHADOZ 
ozonesondes has been made within the inner tropics. For 
the latitude bands of 60-50° N and 50-60° S, representing 
the midlatitudes, trends have been compared between SCIA- 
MACHY and EOS MLS. The trends from SCIAMACHY 
and the comparison instruments are found to be in reason- 
able agreement down to altitudes around 20 km. The shape 
of tropical O 3 trends identified from SCIAMACHY is overall 
well confirmed by all comparison instruments. In the lower- 
most tropical stratosphere below 20 km, SCIAMACHY and 
the comparison instruments see non-negative trends but do 
not agree in detail. The other instruments show higher posi- 
tive trends than SCIAMACHY. 



Appendix A 

QBO and solar cycle variability as seen by SCIAMACHY 

The QBO proxies, the lOhPa and 30hPa Singapore winds, 
are phase-shifted by approximately a quarter of a period. This 
combination can be roughly approximated by the sum of a 
sine and a cosine term with an effective frequency 

Asin(a>f + </>) = A cos (cot) sirup + A sin (cot) cos cp = (Al) 
= a cos (cot) + b sin (cot). 

The phase constant is given by cp — arctan ( |). For t — 0, this 
yields 

Asincp = a, (A2) 

or 


sin (arctan ( | )) 

Determining the QBO amplitude by this approach, a and 
b are given by the fit coefficients of the QBO proxies (cf. 
Eq. 3). 

The inferred QBO amplitude from SCIAMACHY limb O 3 
is shown in Fig. Al as a function of altitude and latitude. At 
equatorial latitudes, a maximum amplitude is obtained at al- 
titudes of 25-28 km and 30-36 km. This pattern of maximum 
amplitude is consistent with the vertical structure of the QBO 
response of tropical O 3 described in Sect. 3.2. Both the spa- 
tial distribution and the strength of the inferred amplitude are 
in close agreement with MIPAS/ENVISAT observations as 
shown by Eckert et al. (2013) who used a similar approach. 
The maximum of SCIAMACHY QBO amplitude of around 
0.6ppmv at around 33 km is larger than MIPAS values by a 
factor of around 1.5. This is consistent with findings by Tegt- 
meier et al. (2013). Comparing QBO signatures measured 
by a variety of instruments, it was shown that the QBO am- 
plitude from SCIAMACHY lies above the multi-instrument 


Fig. Al. Latitude-altitude plot of the amplitude of the QBO re- 
sponse of O 3 as seen by SCIAMACHY: latitude range from 
60° N to 60° S, altitude range from 15 to 50 km. August 2002-April 
2012. Legend on the right: units of ppmv. 
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Fig. A2. Latitude-altitude plot of solar cycle signal in SCIA- 
MACHY limb O 3 : latitude range from 60° N to 60° S, altitude 
range from 15 to 50 km. Significances at the 2cr level are non- 
hatched areas. The solar cycle signal is expressed as the solar max 
minus min change of O 3 (more detailed explanation in the text). 
It is scaled in relative units with respect to the mean value of the 
underlying time series. The period of August 2002-April 2012 is 
considered. Legend on the right: solar max minus min change in 
units of %. 


average, while MIPAS/ENVISAT exhibits a QBO amplitude 
lower than average. 

Some altitude dependence is also discernible in the solar 
cycle signal of O 3 , obtained from the use of the Mg II in- 
dex as solar cycle proxy, shown in Fig. A2. For the Mg II 
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index, a change of 0.01 between solar maximum and min- 
imum is assumed (this is considered to be a realistic esti- 
mate of the difference between the maximum of solar cycle 
24 and the minimum at the end of solar cycle 23). Scaling 
the maximum minus minimum change of the proxy by its 
regression coefficient yields the corresponding solar max mi- 
nus min change in O3. Here, it is expressed in relative units 
(with respect to the mean value of the underlying time se- 
ries). For large parts, the solar cycle signal is significant at 
the 2er level at altitudes between 35 and 45 km and insignif- 
icant below and above. Analysing a merged satellite dataset 
from 1984-2011, the study of Kyrola et al. (2013) inferred a 
largely similar distribution of significances for the solar cy- 
cle signal. The strength of their solar cycle signal being lower 
than ours by a factor up to 2 may be partly explained by the 
difference in the analysed time span and/or scaling the so- 
lar cycle signals in relative units with respect to the constant 
term of their regression (not the mean value of the underlying 
time series). 
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